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^ ■ ABSTRACT 

, Context. AX Per is an eclipsing symbiotic binary. During active phases, deep narrow minima are observed in its light curve, and the 

' ionization structure in the binary changes significantly. From ~2007.5, AX Per entered a new active phase. 

I ' Aims. We aim to derive the ionization structure in the binary and its changes during the recent active phase. 

Methods. We used optical high- and low-resolution spectroscopy and UBVRcIc photometry. We modeled the SED in the optical and 
broad wings of the Ha line profile during the 2007-10 higher level of the AX Per activity. 
pX ' Results. After 10 orbital cycles (~ 18.6 years), we again measured the eclipse of the hot component by its giant companion in the 

light curve. We derived a radius of 27 ±2Rq for the eclipsed object and 115 ±2Rq for the eclipsing cool giant. The new active phase 
was connected with a significant enhancement of the hot star wind. From quiescence to activity, the mass-loss rate increased from 
<~i , ~ 9 X 10"** to ~ 3 X 10"^ Mq yr"' , respectively. The wind causes the emission of the He++ zone, located in the vicinity of the hot star, 

and also is the reason for the fraction of the [O iii] zone at farther distances. Simultaneously, we identified a variable optically thick 
warm (T^ff ~ 6000 K) source that contributes markedly to the composite spectrum. The source was located at the hot star's equator 
and has the form of a flared disk, whose outer rim simulates the warm photosphere. 

Conclusions. The formation of the neutral disk-like zone around the accretor during the active phase was connected with its enhanced 



o 



^ ■ wind. It is probable that this connection represents a common origin of the warm pseudophotospheres that are indicated during the 

^ I ' active phases of symbiotic stars. 

Key words. Accretion, accretion discs - stars: binaries: symbiotic - stars: winds, outflows - stars: individual: AX Per 

1. Introduction To explain the high luminosity of hot components in symbi- 

OO otic binaries, it was suggested that they are powered by stable 

(N| Symbiotic stars are interacting binary systems comprising a hydrogen nuclear- burning on the WD s urface, which requires 

• cool giant and a hot compact star, mostly a white dwarf (WD), ^ certain range of accretion rates (e.g. ' Paczvnski & Zvtkow I 

^ on, typically, a few years' orbit. The WD accretes from the 1978; Puiimoto 1982). During active phases, we observe a 

^ giant's wind, heats up to 1-2x105K, and becomes as lumi- signi ficant coohng of the spect rum produced by the hot star 

^ . nous as 10^ - IQ-^ Lq. It ionizes the circumbinary environ- (e.g. iKenvon & Webbink I flOSi and references therein). The 

. ment and causes nebular emission. As a result, the spectrum event of outbursts could result from an increase in the accre- 

.'^^ of symbiotic stars consists of three basic components of ra- tion rate above that which sustains the stable burning, which 

diation - two stellar and one nebular. If the processes of the jgads to an expansion of the burning envelope to an A- 

^ mass loss, accretion and ionization are in a mutual equilib- p type (pseudo)photosphe re (e.g. iTutukov & Yangelson I [19761 

C3 rium, the symbiotic system releases its energy approximately iPaczvnski & Rudak \ \md). As a result, the pseudophotosphere 

at a constant rate and spectral energy distribution (SED). This radiate at a lower temperature, and consequently shift the 

stage is called the quiescent phase. Once this equilibrium is maximum of its SED from shorter to longer wavelengths, caus- 

disturbed, the symbiotic system changes its radiation signifi- jng ^ brightening in the optical. Based on optical observa- 

cantly, brightens in the optical by a few magnitudes and usually tions, this scenario w as supported by many authors. Recently, 

shows signatures of a mass outflow for a few months to years. LSiviero et al. I (l2509l) applied it to the 2008-09 active phase 

We name fliis stage flie active phase. Many particular aspects of ci Cyg. However, modehng the UV/near-IR continuum of 

of this general view have been originally pointed out by, e.g., symbiotic binaries with a high orbital incUnation during active 

Friediung & Viotti ( 1982), Seaquist, Taylor, & Button ( 1984), phases indicated an increase of both the stellar radiation from a 

iKenyon & Webbinkj il9M), Nussbaumer & Vogel 09871), and warm(l-2xl04K)pseudophotosphereandthenebularradiation 



are more recently discussed in.Corradi et al. . (2003,). within the optical wavelengths. This led to the suggestion that 

there is an edge-on disk around the accretor, whose outer flared 

Tables 2, 3 and 5 are available at http://www.aanda^ rim represents the warm pse udophotosphe re, and the nebula is 

Visiting Astronomer: Asiago Astrophysical Observatory placed above/below the disk (ISkopalll2005h . In this contribution 
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we test this scenario on the recent active phase of the symbiotic 
binary AX Per 

At present, AX Per is known as an ecHpsing symbiotic binarj^ 
(Skonal 1991) comprising a M4.5III giant (Murset & Schmid 
1999.) and a WD on a 680 -d orbit (e.g. Mikolaiewska & Ken von. 



Table 1. Log of spectroscopic observations 



1992t iFekel etal.ll2000h . The last active phase began in 1988, 
when AX Per brightened by ~3 mag in the visual, developed 
a specific phase-dependent modulation in the light curve (LC), 
and showed narrow eclip ses at the position of the inferior con- 
junction of the giant (e.g.'lvison et al. I[l993l) . Transition to qui- 
escence happened during 1995, when the LC profile turned to 
the wave-like orbitall y related variation at a low level of the 
star's brightness ( Skopal etal. I [2001. Fig. 1 here). Recently, 
lMunarietal.1(l2009aF reported on a rapid increase in the AX Per 
brightness by ~1 mag in B during 2009 March. They pointed a 
similarity between this brightening phase and the short-duration 
flare that occuiTed in the AX Per LC about one year be- 
fore it s major 1988-1995 active phase. Recently. iMunari et al. I 
(1201 Ol) reported on a rise in brightness of AX Per during 2010 
November by ~0.7 mag in B, which supported their previous 
suggestion that the 2007-10 active phase was a precursor to an 
oncoming major outburst. From this point of view, understand- 
ing the precursor- type activity can aid us in understanding the 
nature of the Z And-type outbursts better 

Accordingly, in this paper we investigate the structure of the 
hot component in the binary (i.e. the hot star with the warm 
pseudophotosphere and the nebula) that developed during the 
2007-10 brighter phase of AX Per To achieve this aim, we ana- 
lyze our multicolour photometric observations and the low- and 
high-resolution optical spectra. In Sect. 2 we summarize and de- 
scribe our observations and data reduction. Section 3 describes 
our analysis and presents the results. Their discussion and con- 
clusions are found in Sects. 4 and 5, respectively. 



Date 


Julian date 


Res. 


Disp. 


A range 


Telescope 




+2450000 Power (A/pix) 


(A) 




2007/07/31 


4312.876 


10000 




H/3, Ha 


DDO 1.88m 


2007/11/20 


4425.529 




1.81 


3824- 


7575 


CrAO 2.6m 


2007/12/03 


4438.438 




2.30 


3230- 


7770 


Asiago 1.22m 


2008/01/23 


4488.506 


10000 




H/3, Ha 


DDO 1.88m 


2008/07/07 


4655.473 




1.81 


3774- 


7574 


CrAO 2.6m 


2008/08/10 


4689.450 




1.81 


3749- 


7575 


CrAO 2.6m 


2008/09/24 


4734.609 




1.81 


3749- 


7575 


CrAO 2.6m 


2008/10/13 


4753.638 


20 000 




3690- 


7300 


Asiago 1.82m 


2008/10/23 


4763.417 




1.81 


3349- 


7149 


CrAO 2.6m 


2008/11/08 


4778.393 




1.81 


3749- 


7500 


CrAO 2.6m 


2008/12/08 


4809.390 




2.30 


3250- 


7830 


Asiago 1.22m 


2008/12/21 


4822.437 


20 000 




3690- 


7300 


Asiago 1.82m 


2009/01/09 


4841.426 


20 000 




3690- 


7300 


Asiago 1.82m 


2009/03/24 


4915.422 




2.12 


3890- 


7930 


Varese 0.60m 


2009/04/01 


4923.323 




2.30 


4200- 


7200 


Asiago 1.22m 


2009/04/07 


4929.303 




2.12 


3900- 


8120 


Varese 0.60m 


2009/04/08 


4930.301 




2.30 


3340- 


7620 


Asiago 1.22m 


2009/04/09 


4931.262 




1.81 


3600- 


7575 


CrAO 2.6m 


2009/04/14 


4936.281 




2.30 


3560- 


7580 


Asiago 1.22m 


2009/04/14 


4936.338 


17 000 




4030- 


7510 


Varese 0.60m 


2009/05/05 


4957.458 




2.12 


3850- 


8270 


Varese 0.60m 


2009/08/05 


5049.535 




2.30 


3320- 


7660 


Asiago 1.22m 


2009/08/26 


5070.565 




1.81 


3324- 


7493 


CrAO 2.6m 


2009/09/24 


5099.584 




1.81 


3350- 


7500 


CrAO 2.6m 


2009/11/10 


5146.485 


20000 




3690- 


7300 


Asiago 1.82m 


2010/08/07 


5416.368 




2.12 


3910- 


8540 


Varese 0.60m 


2010/08/16 


5425.564 




1.81 


3358- 


7500 


CrAO 2.6m 


2010/08/21 


5430.348 


17 000 




3950- 


8640 


Varese 0.60m 


2010/09/14 


5454.560 




1.81 


3274- 


7574 


CrAO 2.6m 



2. Observations and data reduction 

Our observations of AX Per during its 2007-10 higher level of 
activity were carried out at different observatories. 

Spectroscopic observations were secured with five different 
instruments: (1) the 1.22-m telescope, operated in Asiago by 
the Department of Astronomy of the University of Padova, was 
used with a B&C spectrograph, 300 In/mm grating and ANDOR 
iDus 440A CCD camera, equipped with an LEV 42-lOBU back 
illuminated chip (2048x512 pixels of 13.5 jum size); (2) the 
1.82-m telescope, operated in Asiago by INAF Astronomical 
Observatory of Padova, mounting a REOSC echelle spectro- 
graph equipped with an AIMO E2VCCD47-10 back-illuminated 
CCD detector (1 100x1 100 pixels of Ujjm size); (3) the 0.60- 
m telescope of the Schiaparelli Observatory in Varese, mount- 
ing a multi-mode spectrograph, able to provide both single dis- 
persion, low-resolution spectra and Echelle spectra. The detec- 
tor was an SBIG ST-IOXME camera (2194x51472 pixels of 
6.8 i-im size); (4) the 1.88-m telescope of the David Dunlap 
Observatory (DDO), University of Toronto (DDO), equipped 
with a single-dispersion spectrograph equipped with a Jobin 
Yovon Horiba CCD detector (2048x512 pixels of 13.5 yum size); 
and finally (5) the 2.6-m Shajn telescope, operated by the 
Crimean Astrophysical Observatory (CrAO) in Nauchny, mount- 
ing a SPEM spectrograph in the Nasmith focus. The detector was 
a SPEC-10 CCD camera (1340x100 pixel). Table 1 provides 
a journal of the spectroscopic observations. At each telescope 
the observations were carried out as multiple exposures to avoid 
saturation of the strongest emission lines. Various spectropho- 
tometric standards were observed each night to flux into abso- 



lute units the spectra of AX Per. The accuracy of the flux scale 
was checked against the photometric observations by integration 
over the fluxed spectra of the UBVRcIc photometric pass-bands. 
Correction for flat-, bias- and dark frames was caiTied out in a 
standard way with IRAF, as well subtraction of the background 
sky and scattered light. 

Multicolour BVRqIc CCD photometry was obtained with 
various telescopes operated by the ANS (Asiago Novae and 
Symbiotic Stars) Collaboration. Treatments for flat-, bias- and 
dark frames was caiTied out in a standard way. Photometric cali- 
bration and correction for color equations wa s carried out against 
the p hotometric sequence calibrated by iHenden & Munaril 
(I2006h around AX Per. The resulting data are presented in 
Table 2. Their uncertainties are the overall budget errors (which 
includes the Poissonian component and the transformation to the 
standard Landolt system). They are of a few xO.Ol - 0.001 mag. 

In addition, classical photoelectric U BV measurements were 
performed by a single-channel photometer mounted in the 
Cassegrain focus of 0.6-m reflector at the Skalnate Pleso ob- 
servatory. Internal uncertainties of these one-day-mean mea- 
surements are of a few xO.Ol mag. These data are presented in 
Table 3. 

Arbitrary flux units of the high-resolution spectra were 
converted to absolute fluxes with the aid of the simultaneous 
UBVRcIc photometry. To determine flux-points of the true con- 
tinuum from the measured mag nitudes, we ca lculated coiTec- 
tions for emission lines. Am/ (see lSkoparil2007h . using our low- 
resolution spectra (Table 4). Then we interpolated these values 
to dates of the high-resolution spectra. S pectroscopic observa- 
tions were dereddened with Eb-v - 0.27 dKenvon & Webbiiikl 
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Fig. 1. Top panel displays the historical LC of AX Per from ~1890 (see Fig. 1 of lSkopal et al. 1120011 and references therein). We 
used visual magnitude estimates from the AAVSO International Database and those gathered by members of the AFOEV, which are 
available on CDS. Visual data were smoothed within 20-day bins. Bottom panel shows the UBV LCs of AX Per from its last 1988 
major outburst. Minima timing corresponds to the ephemeris given by Eq. (1). 



|1984|) an d the resulting parameters were scaled to a distance of 
1.73 kpc (ISkot^al 1120001) . 



3. Analysis and results 

3. 1 . Photometric evolution 

3.1 .1 . Indications of a new active pliase 

According to the LC evolution from 1988 (Fig. 1), the increase in 
the star's brightness by AU ~Q.5 mag from 2007 July signals that 
AX Per entered a new active phase. An additional rapid bright- 
ening by AB ~0.8 mag and a bluer index B -V ~ 0.8, observed 
during 2009 March (Munari et al. 2009a, and Fig. 2 here), sup- 
ported the new activity of AX Per A significant change of the 
broad minima, observed during quiescent phases at/around the 
giant inferior conjunction, into a naiTower and deeper eclipse, 
indicates the active phase of a symbiotic binary as well (e.g. 
^elyakina 1979, 1991, and Figs. 1 and 2 here). During the 
2007 conjunction we measured a relatively small, V-type min- 
imum, which position was shifted from the giant conjunction 



by ~ -0.025 Poib- During the following conjunction the eclipse 
was deeper and broader, nearly rectangular in profile and placed 
at the inferior conjunction of the giant. The light observed dur- 
ing the totality became redder with B - V ~ 1.4 (see Fig. 2). 
Finally, we measured a light wave in the LC with a period of 
~ 0.5 f orb and minima located around orbital phases 0.2 and 0.7 
(Fig. 2, left). Its presence was transient, being connected only 
with the 2007-10 higher level of the star's activity (see Fig. 1). 
It can be seen in all filters, but especially in the Iq band is well 
pronounced also after the 2009 eclipse (see Fig. 2). This feature 
was observed in a symbiotic system for the first time. Its un- 
derstanding requires a detailed analysis, which is, however, out 
of the scope of this paper. Therefore we postpone it to a future 
work. 



3.1.2. Eclipses 

The minimum (eclipse), observed during the 2009 brightening, 
was very similar in profile to those observed during the major 
1988-91 active phase. Its well pronounced profile closely cov- 
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Fig. 2. Left: Recent UEVRqIc LCs of AX Per covering its 2007-10 active phase. Thin and thick arrows denote the dates of our 
low- and high-resolution spectra, respectively (Table 1). Vertical dotted lines mark the time of eclipses according to Eq. (1). Right: 
Minima profiles selected from LCs in Fig. 1, and folded with the ephemeris (1). Vertical dashed and solid line mark the light 
minima during quiescent and active phases, respectively. Denotation in keys, for example, "Q: U4-2.2", means quiescent phase, U 
filter, epoch 4 and a shift by -2.2 mag. 



ered by measurements, mainly in the B- and V-bands, allowed us 
to determine its contact times as Ti - JD 2454998.50 + 0.79, 
T2 = JD 2455 012.50 + 0.53, T3 = JD 2455 065.40 + 0.85 
and T4 = JD 2455 085.33 + 0.50, where uncertainties include 
only those given by the data coverage, assuming a linear de- 
pendence of the star's brightness during the ingress to and as- 
cent from the totality. We neglected other possible sources of 
errors. These times and their uncertainties indicate an asym- 
metric profile with T2 - Ti = 14 + 1.3 d and T4 - Tj, - 
20 + 1.3 d, which implies an increase of the linear size of the 
eclipsed object during the totality. Therefore, we determined the 
middle of the eclipse as the average of all contact times, i.e. 
7Dec1.(2009) ^JD2 455 040.43 + 0.67. Using well-defined tim- 
ings of other eclipses, 7Dec1 (1989) = 7Z) 2 447 55 1 .7 + 1.0 and 
7Dec1.(1990) = 7D 2448 231.6 + 0.6 (ISkoDallll99lh . gives their 
ephemeris as 

JDeci. = 2 447 551.26(+0.3)-h680.83(+0.11)x£. (1) 



Our photometric ephemeris agrees within uncertainties with that 
given by the solution of the spectroscopic orbit, Tsp conj. = 
2 447 553.3(+5.6) -1- 682.1(±1.4 ) x E obtained from infrared 
radial velocities by iFekel et aTl (l2000l) . Finally, assuming the 
orbital inclination ; - 90° (see a summary in 'Skopa l et al. I 
2001) and the separation between the binary components, a = 
364/^0 ( from the mass ratio and the mass func tion published 
by Mikolaiewska & Kenvon 1 1 1 9921: iFekel et al. 1 ,2000) . the ra- 
dius of the giant, Rg = 1 15 ±2 Rq and that of the eclipsed object, 
Rf. - 27 + 2Rq. This result suggests an expansion of the giant 
radius by ~13^0, from 102 + 3R(7), given by the contact times 
of the 1990 eclipse (ISkopallll994 . The asymmetry of the 2009 
eclipse profile and the radius of Tie = 36 + 3 Rq derived from the 
1990 eclipse profile suggest that the eclipsed object is subject to 
variation (see also Table 8 below). 
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-250 250 -250 250 -250 250 -250 250 -250 250 

Heliocentric radial velocity [km s''] 

Fig. 3. Evolution in the selected line profiles along the 2007-10 active phase as observed in our high-resolution spectra. Fluxes are 
in units of 10"'^ erg cm"^ s"' A"' . Systemic velocity of -1 17.44 km s"' was subtracted. Small numbers at the right side of panels 
represent a shift of profiles with respect to the level of the local continuum. The spectra are dereddened. 



Table 4. Effect of emission lines on the UBVRc magnitudes 



Date 


KU, 


AS, 


AK, 


KRcj 


2007/11/20 




-0.26 


-0.07 


-0.20 


2007/12/03 


-0.14 


-0.24 


-0.05 


-0.15 


2008/07/07 


-0.10 


-0.30 


-0.06 


-0.20 


2008/08/10 


-0.11 


-0.33 


-0.07 


-0.22 


2008/09/24 


-0.12 


-0.34 


-0.07 


-0.24 


2008/10/23 


-0.20 


-0.35 


-0.08 


-0.25 


2008/11/08 




-0.34 


-0.07 


-0.24 


2008/12/08 


-0.18 


-0.32 


-0.06 


-0.21 


2009/03/24 




-0.19 


-0.05 


-0.13 


2009/04/01 




-0.15 


-0.06 


-0.12 


2009/04/07 




-0.17 


-0.05 


-0.12 


2009/04/08 




-0.20 


-0.05 


-0.12 


2009/04/09 




-0.21 


-0.06 


-0.14 


2009/04/14 




-0.19 


-0.06 


-0.12 


2009/05/05 




-0.14 


-0.04 


-0.10 


2009/08/05 


-0.15 


-0.23 


-0.06 


-0.10 


2009/08/26 


-0.10 


-0.21 


-0.07 


-0.13 


2009/09/24 




-0.15 


-0.06 


-0.13 


2010/08/07 




-0.30 


-0.07 


-0.22 


2010/08/16 


-0.18 


-0.36 


-0.08 


-0.27 


2010/09/14 


-0.22 


-0.39 


-0.10 


-0.31 



3.2. Spectroscopic evolution 
3.2.1. The line spectrum 

During the investigated period of the 2007-10 activity, the op- 
tical spectrum of AX Per was dominated by Hi, He i and He ii 



emission lines, although some faint emissions, produced by 
permitted transitions in ionized metals, as e.g. Niii4641 and 
Cm 4647 A, were also present. The spectrum of the forbid- 
den lines was characterized mainly by nebular lines of [Oiii], 
[Neiii] and two [Fevii] lines at 5721 and 6087 A. The Raman- 
scattered O VI 1032 line at 6825 A was not present at all. Figure 3 
shows evolution in the line profiles of Hen 4686 A, li/3, Ha, 
[O III] 5007 A and [Fe vii] 6087 A, as observed on our high- 
resolution spectra. Figure 4 demonstrates variations in the line 
fluxes along the active period, which we observed in all our 
spectra. Their quantities are introduced in Table 5 (available on- 
line). We also measured radial velocities (RVs) of main emission 
peaks in the line profiles of highly ionized elements (Table 6). 
Figure 5 displays their values compared to the radial velocity 
curves derived by other authors. 

Hydrogen lines: Profiles of hydrogen B aimer lines were sig- 
nificantly affected by a blue-shifted absorption component. Its 
presence can be recognized around the whole orbital cycle. From 
the orbital phase (f = 0.85 to ip - 1.19, it created a double- 
peaked profile (Fig. 3). In contrast, during the previous quiescent 
phase, the Ha line observed at similar orbital phases {tp - 0.81 
and ip = 0.16, as on 1998 Jan. 9 and Sept. 5, respectively, 
see Fig. 5 in Skopal et al. 2001), did not display the double- 
peaked profile. In addition, both the total flux and the extension 
of the broad wings of mainly the Ha profile were significantly 
higher than those measured during quiescent phase (Fig. 6 and 
Sect. 3.2.2 below). It is important to note that during the ma- 
jor 1988-91 active phase of AX Per a strong absorption cut the 
broad emission profi le even at position s with the hot star in front 
(see Figs. 5 and 8 of llvison et al. 1fT99"3ft . This indicates the pres- 
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ence of neutral material along the line of sight around the whole 
orbit, which attenuates line photons via the b - b transitions. 
The strong difference in the line profiles measured during ac- 
tive and quiescent phases suggests that the density of the absorb- 
ing material enhances during active phases at the orbital plane. 
Finally, during the eclipse the total fluxes of Ha and lines and 
those of their wings decreased by a factor of > 3 with respect 
to their maximum values (Fig. 4, Table 5). The s ame behavior 
was o bserved during the 1994 eclipse (see Fig. 6 of lSkopal et al. I 
12001 ). Hence, during the active phases, most of the ionized re- 
gion around the hot star is confined in a region that has about the 
size of the red giant's radius. 

Heii4686 A line: The Heii4686A emission line nearly dis- 
appeared during the echpse. This locates the region of its origin 
in the hot star vicinity and constrains its maximum linear size 
to the diameter of the eclipsing giant. The measured flux in the 
He II 4686 line, f 4686, allows us to estimate the mean electron 
concentration in the He^^ zone. Its luminosity can be expressed 
as 

4;r^^^f4686 = Q'4686(7'e)nen(He++)/lV4686V4686, (2) 

where a4(,^(,(Tg) is the effective recombination coefficient for the 
given transition, and n(He^^) is the mean concentration of 
electrons and He^^ ions, respectively, and V4686 is the volume of 
the He^^ zone. We also assume that all He atoms are fully ion- 
ized within the zone, i.e. the abundance A(He^^) - A(He), and 
thus n(He^^) - A(He )np, because the concentration of hydrogen 
atoms is equal to that of protons, Wp, in the ionized zone. Then 
the electron concentration 

ne = (l+2A(He++))«p, (3) 

because each He atom produces two electrons within the He^^ 
region. For A(He^^) = 0.1 the ratio np/ng ~ 0.83, and Eq. (2) 
can be expressed as 

4nd^F46S6 = (^^4^) t^4686(re)A(He )«^0.83/jV4686, (4) 

where we approximated the V4686 volume by a sphere with 
the radius of the eclipsing giant, R„, and e is the filling fac- 
tor Our measured fluxes, ^4586 ~ 1.56 x 10"" and ~ 7 x 
10 erg cm"^ s produced by the He^^ zone before and after 
the eclipse, respectively (Fig 4, Table 5), a4686(3 000 K)/iV4686 
= 4.2 X 10-2^ ergcm^s ' dHummer & Storey 1 1 1987). = 
USRq (Sect. 3.1.2.) and e = 1, yield a mean = 8.6 x 10'' 
and 5.7 X 10^ cm"^ during the 2009 brightening, before and af- 
ter the eclipse, respectively. These quantities represent a mini- 
mum rig, because an optically thin case is assumed. The size of 
the He^^ zone can also be estimated from spectroscopic obser- 
vations made close to the Tt, time (spectrum 2009/08/26) and 
just after the T4 time (2009/09/24), when the He 11 4686 line also 
aroused from the eclipse (Fig. 4). This time interval (~30 days) 
corresponds to /?4686 ~ 50 Rq, which yields ~ 2.0xl0'"cm"^ 
(volume of the eclipsed object was not subtracted). 

Profile of the He 11 4686 A line was not symmetrical with re- 
spect to the reference wavelength. The RV of its peak was placed 
mostly at the antiphase to the RV curve of the giant (Fig. 5), 
which supports the connection of the He^^ zone to the hot star 
surroundings. Its larger amplitude indicates an extension of the 
He^^ zone to farther distances from the hot star. However, the 
position of the main core of the profile, measured around the half 
of its maximum, was shift ed blueward relative to the sy stemic 
velocity (see also Fig. 5 of iMikolaiewska & Kenvonll992) . The 



Table 6. Radial velocities in kms ' of the main emission peak 
of highly ionized elements. The systemic velocity was sub- 
tracted. 



Date 


He II ^4686 


[Fevii] ^6087 


[Oiii] .15007 


2007/07/31 


10.5 






2008/01/23 


2.0 






2008/10/13 


6.9 


10.3 


27.1 


2008/12/21 


20.3 


24.8 


33.4 


2009/01/09 


22.7 


25.8 


31.7 


2009/04/14 


23.8 




-16.2 


2009/11/10 


-18.0 




-14.1 


2010/08/21 


-0.6 


-13.4 


29.4 



violet shift of the He 11 4686 A emission line was measured also 
in the spectra of other symbiotics. For example, Tamura et aO 
(12003b i neasured it fo r Z And and Munari et al. I ( 12009b() for 
AG Dra. lSkopal et al. I (I2OO61) interpreted this efi'ect by the pres- 
ence of a disk around the hot star, which blocks a part of the 
radiation from the wind at its backside to the direction of the ob- 
server. Finally, the profile seems to be affected by a variable ab- 
sorption from its violet side, similarly as for the hydrogen lines 
(Fig. 3), which makes it difficult to interpret the He 11 4686 A line 
profile more accurately. 

Forbidden lines: In contrast to the H and He fines, the neb- 
ular A^i and N2 ([Oiii] 5007 and 4959 A) fines were not subject 
to the eclipse. Fluxes of the stronger 5007 A line were as low 
as ^ 1 X 10 erg cm"^ s"' before 2009 March, when they in- 
creased to ~ 4.5 X 10"'^ erg cm"^ s"\ following the increase 
in the optical continuum. However, during the eclipse the line 
flux persisted at the same level. In 2010 it returned to the pre- 
brightening values (Fig. 4). The ratio R - F(Ni + A^2)/-F'4363i 
which is a well-known probe of We and in planetary nebu- 
lae, was extremely low. The observed values, 0.5 - 2, imply a 
high-density [Oiii] nebula withne([0 111I) ~ 10^ - 10^ cm~^ and 
re([Oiii]) < 20000K (e.g. Fig. 15.1 in lGurzadvan]ll997h . A 
detailed application of the method to AX Per can be found in 
ISkopal et alTI (12001 ). who derived the upper limit of ne([Oiii]) 
as 7 X 10^ cm"^ for R = 1.7 - 4.3. It is interesting that during 
the 1994 eclipse fluxes in the A^i and N2 lines faded by a factor 
of <2 with respect to their out-of-eclipse values (see their Fig. 6 
and Table 4), in contrast to the evolution during the 2009 eclipse. 
We discuss this effect in Sect. 4.1.1. 

The highly ionized [Fe vii] 6087 A line was present in the 
spectrum only during a lower activity level. From the opti- 
cal brightening in 2009 March to our observation on 2009 
November 10th, it entirely disappeared from the spectrum 
(Table 5). Its profile was asymmetrical with respect to the ref- 
erence wavelength, being broadened more to the violet side 
of the line (Fig. 3). As a result, the RV of the intermediate 
point at half of the line intensity was negative relative to the 
systemic velocity (-18.1, -9.0, -7.3 and -13.4km s"' as on 
2008/10/13, 2008/12/21, 2009/01/09 and 2010/08/21, respec- 
tively). However, the peak of the emission line followed the RV 
curve as derived by lijima ( 1988); (see Fig. 5 here), which im- 
plies that a part of the [Fe vii] zone has to be connected with 
the hot star, but it is located farther from it, behind the bi- 
nary. The [Fe vn] transitions are excited by colli sions with ther- 
mal electrons (iNussbaumer & Osterbrockl 19701) . The enhanced 
ionized wind during the brightening produces a larger amount 
of free electrons, while the [Fe vii] 6087 A emission line disap- 
pears. This anti-correlation can be understood if the [Fe vn] re- 
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gion is placed instead at the orbital plane, where the material 
is in the ionization shadow during the phase of a warm disk 
creation (see Sect. 4). During quiescence, when the disk is not 
present, the [Fevn] line can be created withi n the coUisional 
zone of winds from the binary components ( Wallerstein et al. I 
19841). According to its shaping in the binary (see Fig. 5 of 



example, M„ = 3 x 10 ^M0yr ' and - 2000 km s \ yield 



Bisikalo et al. 112006'). a fraction of the [Fe vii] emission could 



arise behind the binary from the hot star side, and thus sat- 
isfy the RV of its orbital motion with a larger amplitude (see 
Fig. 5). However, more high-resolution observations of the 
[Fe vii] 6087 A line for eclipsing symbiotic binaries are needed 
to understand its origin in symbiotic systems better. 

3.2.2. Mass loss through the hot star wind 

The hot components in symbiotic binaries can lose their mass 
in the form of wind (e.g. Vogel 1993t Nussbaumer et al. 11199^1 
[Eriksson et al. 1 120041: iKennv & Tavlor I 2005h. wh ich enhances 
considerably during active phases (ISkopalll2006l). To estimate 
its rate, we used the Ha method as proposed by Skopal (2006), 
which assumes that the broad Ha wings originate in the ionized 
wind from the hot star. The principle of the method is a compar- 
ison of the observed and synthetic profile of the broad wings. 

The model assumes a spherically symmetric wind originat- 
ing in the central star, which is covered by an optically thick disk 
in the direction of the observer. The disk is characterized with the 
radius Ru and height Ho, and it is seen edge-on in the model (see 
Fig. 9). Density, and thus the emissivity of the wind at a given 
distance r from the hot star centre, is determined by the mass- 
loss rate and the veloc ity of the wind. Th e velocity distribution 
is approximated with (iCastor et al. Ill975h 



(5) 



where the origin of the wind, R^, and /3 are model parame- 
ters, while the terminal velocity, Vco, is given by the extension 
of the wings. During the quiescent phase, parameters Ru and 
//d can be estimated from the effective radius of the hot star, 
Rf ~ 0.1 Rq (Table 3 of Skopal 2005). Assuming the ratio 
Hd/Rd = 0.1 yields Ru = 0.32 ^0 and Hq = 0.032 Rq. During 
the investigated active phase, we estimated the disk radius from 
the first two contact times of the 2009 eclipse (Sect. 3.1.2.) as 

~ 20 Rq. Assuming a disk with Hu/Rd = 0.3, yields Hu - 
6 Rq. We recall here that the results do not critically depend 
on these parameters. The mass-loss rate is mainly determined by 
the luminos ity of the broad wings and their terminal velocity (see 
Eq. (14) of lSkopalll2006h . A comparison between the modeled 
and observed profiles is shown in Fig. 6. The synthetic profiles 
match the observed wings well for \RV\ ^ 200 - 250km s"'. 
There is a significant difference between the extension and lumi- 
nosity of the Ha wings observed during the quiescent and active 
phase (Table 7, Fig. 6). Broad wings from the 1998 quiescence 
corresponded to the mass-loss rate of ^ 1 x 10"^ Mq yr"', while 
during the activity it increased to 2 - 3 x 10"^ Mq yr"'. 

In our model the mass-loss rate determines the particle den- 
sity at a given distance r from the origin of the wind according 
to the continuity equation as 



n^(r) = Mw/47rr ijmiiv(r). 



(6) 



where the wind velocity, v(r), is given by Eq. (5), /v is the mean 
molecular weight and mn is the mass of the hydrogen atom. 
Within the He^^ zone, i.e. from the hot star to rHeii ~ 50 Rq 
(Sect. 3.2.1), the particle densities n„{r) excellently agree with 
those derived independently of the measured F4686 fluxes. For 



hMORq) ~ 2.8 X 10'" (i.e. = nw/0.83 = 3.3 x 10"' cm 
and n„(5O/?0) ~ 3.2xlO''cm-3 (n^ = 3.9 x lO'^cm"^). This sug- 
gests that particles of the hot star wind give rise to the observed 
emission from the He^^ zone. 



3.2.3. Modeling the continuum spectrum 

The observed continuum spectrum of symbiotic stars, F(A), can 
formally be expressed as a superposition of its three basic com- 
ponents (see Sect. 1), i.e. 



F{A) = F^(A) + F,M) + Fi,(A), 



(7) 



where Fg{A), Fn(A) and Fh(A) represent radiative contributions 
from the cool giant, nebula and a hot stellar source, respectively. 
Basic features of these components of radiation are well recog- 
nizable in our low-resolution spectra exposed from short wave- 
lengths, A ^ 3300 A. Most distinctive are molecular passbands 
arising in the red giant atmosphere, while a pronounced Balmer 
jump in emission directly indicates a strong nebular continuum. 
The source of ionizing photons (i.e. the hot star) is indicated 
indirectly by the presence of strong recombination lines super- 
posed on the continuum. Its direct indication is possible only 
within the super-soft X-ray and fa r-UV; in the optical its contri- 
bution is negligible (e.g. Fig. 2 of ISkopal et al. 1120091) . Another 
significantly cooler stellar source of radiation develops during 
active phases. It is connected with the hot star, and its evidence 
is given directly by eclipses in the optical LCs (e.g. Belvakina] 
199li) and/or by model SED at any other orbital 



119791: ISkopal 



phase (ISkopalboOSi ). Hereafter we call it a warm steUar compo 



nent (WSC). 

The aim of this section is to disentangle the composite spec- 
trum, observed at different stages of AX Per activity, into its in- 
dividual components of radiation, i.e to determine their physical 
parameters. For the Fg{A) component we us ed photospheric syn - 
thetic spectra of M-giant stars published bv lFluks et al. I (Il994l) . 
These models were calculated in the spectral range 350-900 nm 
for 11 spectral types (ST), from MO to MIO. In our modeling 
we determine the ST and its subclass by a linear interpolation 
between the neighboring best-fitting STs. Because the radiation 
from the giant can vary, the scaling factor of the synthetic spectra 
represents another parameter in the model SED. 

We approximated the spectrum produced by the nebula by 



F„(A) ^k„x sAT,), 



(8) 



where the scaling factor [cm^^] represents the emission 
measure (EM) of the nebula scaled with the distance d, i.e., 
EM = 4nd^ X k„. e^(re) is the volume emission coefficient 
[ergcm^ s"^ A '], which depends on the e lectron temperature, 
Te, and is a function of the wavelength (e.g. lBrown & Mathews I 
1970). By this approach we can register only those photons pro- 
duced by the optically thin part of the nebular medium, and 
thus estimate the lower limit of EM. For the sake of simplic- 
ity, we calculated the e^(re) coefficient for the hydrogen plasma 
only, including contributions from recombinations and thermal 
bremsstrahlung. Some addition al arg uments for this approxima- 
tion are given in ISkopal et aP (120091) . In addition, we assumed 
that and thus sa{Ts) are constant throughout the nebula. 

The presence of F-G type features in the spectra before and 
after the 2009 eclipse suggests that the temperature of the stellar 
source is as low as » 5000 - 7000 K. As the most obvious F- 
G features we identified the broad Ca 1/14227 and CaiiK/13934 
line (Fig. 7). Also the absorption at /}38 10, pronounced mainly in 
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Table 7. Parameters of synthetic models of the broad Ha wings {R„, /3, Voo), their luminosities observed for \RV\ > 200km s ', 
L„(200), and corresponding mass-loss rates, M„. 



Date 


Rw 






L„(200) 


log(M„) 


dd/mm/yyyy 






[km s-'] 




[Moyr"'] 


Quiescent phase: Rd = 


Q32Rq,Ho = 0.032 


O 


1998/01/09 


0.027 


1.75 


1500 


0.90 


-7.04 


1998/09/05 


0.030 


1.75 


1600 


0.88 


-7.02 




Active phase: Ro 


= 20Rq,Hd 


= 6Ro 




2007/07/31 


5.20 


1.75 


2500 


2.97 


-5.52 


2008/01/23 


5.27 


1.74 


2500 


2.13 


-5.59 


2008/10/13 


5.40 


1.75 


1500 


0.83* 


-5.84 


2008/12/21 


5.45 


1.72 


2000 


1.60 


-5.70 


2009/01/09 


5.40 


1.75 


2000 


1.40 


-5.72 


2009/04/14 


5.25 


1.73 


2500 


3.17 


-5.49 


2009/11/10 


5.20 


1.75 


2100 


1.42* 


-5.67 


2010/08/21 


0.94* 


1.80 


2000 


2.25* 


-5.96 



* \RV\ = 250kms-i,*i?D = 3.5Rq,Hd = 1.05 



the 2009/04/09 spectrum, is clearly evident i n the F and G spec- 
tral tem plates showed on Figs. 218 to 222 of iMunari & Zwitterl 
(l200l . Therefore, we compared the continuum of the hot stel- 
lar source in Eq. (7) to synthetic s pectra, calculated for 5 000 - 
10 000 K using the Kurucz's codes dMunari et al. Il2005h . So, the 
third term of Eq. (7) can be written as 



(9) 



where T^AiTeff) is the WSC of radiation. Its coiTesponding elfec- 
tive temperature, T^ff, and the scaling factor, fl^*^;-,^ - R'^^^^^/d, 
represent free parameters in modeling the SED. Other atmo- 
spheric parameters of synthetic spectra were fixed (log(g) - 
2.5, [M/H] = 0,[Q'/Fe] = 0,Vi-ot = 20kms-i), and their reso- 
lution was accommodated to that of our low-resolution spectra. 

After defining the individual components of radiation in 
Eq. (7), we prepared a grid of models for reasonable ranges of 
the fitting parameters (ST, its subclass and scaling for the gi- 
ant; Te and for the nebular continuum; T^ff and 0^^^.^ for the 
WSC), and selected that which corresponded to a minimum of 
the function 



X 



i=l 



AF°^HAi) 



(10) 



where F^'^'iA-,) are fluxes of the observed continuum, is their 
number (1500 - 1800), AF°^''(Ai) are their errors, and F{Ai) are 
theoretical fluxes given by Eq. (7). Input flux-points were se- 
lected from the observed spectra by omitting emission lines and 
a spectral region from 3645 to ~3740 A, where hydrogen lines of 
the high members of Balmer series were blended, which did not 
allow us to identify the true continuum. For the flux uncertainties 
we adopted typical values of 10%. 

Our resulting models of the AX Per composite continuum 
along its 2007-10 active phase are shown in Fig. 8 and the 
coiTesponding parameters are given in Table 8. The presence 
of the WSC in the 3200-7500 A spectrum is constrained by a 
low Balmer jump and high values of the observed fluxes in the 
blue part of the spectrum. To fit the small Balmer discontinuity 
and to fill in the large difference between the contribution from 
the giant and the observed spectrum for A 5 5000 A, an addi- 
tional radiative component with a fairly flat energy distribution 
throughout the optical was required. This component satisfies 



the radiation from a stellar source with temperature of 5000- 
6000 K well, and can accordingly be associated with the WSC. 
The WSC and the nebular emission were strongest during the 
2009 brightening (Table 8). The following eclipse interrupted the 
brightening by removing the WSC entirely and the nebular com- 
ponent partially from the spectrum (EM decreased by a factor of 
^3). Immediately after the eclipse, the WSC arose again in the 
spectrum together with a stronger nebular emission (Fig. 8 and 
Table 8). We were able to fit our last spectra from 2010 August 
and September without the WSC, which signaled the return of 
AX Per to quiescent phase. For comparison, we inserted the last 
spectrum into the U V-IR SED fro m quiescence (adapted accord- 
ing to Fig. 21 of Skopal 1 120051) to demonstrate the negligible 
contribution of the hot star in the optical (Fig. 8, bottom). 

For comparison, we also derived T^, EM and the intrinsic V- 
magnitude of the giant, Vg, from the observed UBV magnitudes 
by using the method of Carikova & Skopal ( 2010). For the 2009 
eclipse we used B = 13.01 and V = 11.59, taken simultane- 
ously with the low-resolution spectrum from 2009/08/05. The U 
magnitude was interpolated to 12.93 at this date. After simul- 
taneous coiTections for emission lines (Table 4), we determined 
Te ~ 37000 K, EM = 2.7 + 0.2 x lO^** (ii/1.73 kpc)^ cm^^ and 
Vg = 10.98 + 0.02. Uncertainties of the resulting parameters 
reflect errors of photometric measurements. The Vg flux-point 
matches the contribution from the giant well (Fig. 8), the EM 
is also well comparable with that derived from the model SED 
(Table 8), but Tg suffers from a large uncertainty caused by a 
weak dependence of the nebular continuum profile on within 
the U BV region only. 



4. Discussion 

From 2007 on, AX Per entered a new active phase (Sect. 3.1.1). 
This was also supported by a coiTelation between the visual 
brightness, EM and the WSC luminosity (see Table 8), as ob- 
served for the hot components during active phases in other sym- 
biotic stars (see Table 4 of Skopal I I2005I) . The 2007-10 active 
phase was connected with a significant change in the ionization 
structure in the binary during the transition from the preceding 
quiescent phase. 

(i) In the continuum, this was indicated by a distinctive 
change in the shape of the minimum observed at/around the infe- 
rior conjunction of the giant. During quiescence, the broad min- 
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1 ' r 

AX Per 



2007/12/03 



s 

X 




2008/12/08 
Minimum 




-1 

2 

1 h 



2009/04/09 
Brightening 



2009/08/05 
Eclipse 



2009/09/24 



Vp magnitude Q 






3 

E 

o 

^ 



Quiescent phase 2010/09/14 



Observation — 
IJHKLM • 

Fg(^) 

Fh(^) — 
model SED — 




lUE: 1984/10/23 

I 1 I 



3.8 4 
log(Wavelengtli) [A] 

Fig. 8. Our low-resolution dereddened spectra (violet lines) and their models (heavy black lines) at selected dates during the 2007- 
10 active phase of AX Per The model SEDs and their components of radiation here represent a graphic form of Eq. (7) with the 
same denotation in keys (in the bottom panel). Fluxes are in units ergcm^^s ' A Corresponding parameters are collected in 
Table 8. Modeling is described in Sect. 3.2.3. 
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Table 8. Parameters of the SED-fitting analysis (see Sect. 3.2.3, Fig. 8). 



Date 


Stage/B-mag 




Giant 






Warm stellar object 




Nebula 


T 


yyyy/mm/dd 




ST 


T'eff/K* 




Teff/K 


^warm/% 


^warm / ^0 


TJK 


EM/IO^' cm"^ 


2007/12/03 


A / 12.5 


M5.8 


3316 


1430 


5250 


6.2 


26 


34000 


3.5 


0.75 


2008/10/23 


A / 12.4 


M5.8 


3316 


1430 


5250 


8.2 


45 


23 000 


3.8 


1.17 


2008/12/08 


A / 12.7 


M5.8 


3316 


1430 








32000 


4.8 


0.87 


2009/04/09 


A/ 11.6 


M5.7 


3326 


1450 


5500 


11 


98 


25 000 


8.2 


0.80 


2009/08/05 


A / 13.0* 


M5.8 


3316 


1430 








33 000 


2.6 


0.78 


2009/08/26 


A / 12.8* 


M5.6 


3336 


1470 








28 000 


1.7 


1.06 


2009/09/24 


A / 12.0 


M5.6 


3336 


1470 


6250 


7.0 


68 


32000 


4.0 


0.55 


2010/08/16 


Q/ 12.4 


M5.3 


3366 


1520 








27 000 


4.0 


1.30 


2010/09/14 


Q/12.5 


M5.3 


3366 


1520 








27 000 


3.5 


1.67 



* according to the lFluks et al. I ( fT994h calibration, * foriJg = 115i?0' * eclipse 



ima in the LC (e.g. E = 4 to 9 in Fig. 1) are caused by the orbital 
motion of extended ionized wind from the giant, which is par- 
tially optically thick dSkopal 11200 lb . During the active phase, the 
broad minima were replaced by naiTow ones (eclipses), which 
suggests that a large part of the optical light produced by the hot 
component was concentrated at/around the hot star 

(ii) In the line spectrum, we indicated a significant increase 
in the mass-loss rate via the hot star wind with respect to the qui- 
escent phase (Sect. 3.2.2, Fig. 6, Table 7). For example, during 
the 2009 brightening, we measured the highest luminosity and 
extension of the Wa wings, which corresponded to the mass-loss 
rateof ~ 3 x lO^^Moyr"'. 



4. 1 . Evolution of the ionized zone 



Throughout the whole 2007- 10 period we observed a strong neb- 
ular component of radiation in the spectrum (Fig. 8, Table 8). 
A geometrical change of the main source of the nebular radia- 
tion was indicated by the change in the position and profile of 
eclipses during the 2007 and 2009 inferior conjunction of the 
giant. The 2007 minimum was shifted from the conjunction by 
0.025f oib (Fig. 2, right). This implies that an enhanced emit- 
ting region had to be located in fro nt of the hot star, prec eding 
its orbital motion. As discussed by iBisikalo et al. I (l2006l) . this 
higher density region can develop about 100-150 days after an 
increase of the hot star wind as a result of its collision with the 
giant's wind and the orbital motion (see their Figs. 4 and 5). This 
region will be ionized by the neighboring hot star, and thus pro- 
duce a surplus of nebular emission. As a result, the light mini- 
mum will precede the inferior conjunction of the giant. However, 
after the following orbital cycle, the minimum was located ex- 
actly at the giant's conjunction (Fig. 2). Changes in the position 
of the light minima during different stages of activity were in- 
dicated also for other eclipsing symbiotic binaries (see .Skopal i 
[T998h . 

During the 2009.6 eclipse the WSC of radiation disappeared 
entirely, and the nebular component faded by a factor of ~3 with 
respect to the out-of-eclipse values (Table 8, Fig. 8). This im- 
plies that a strong source of the nebular continuum had to be 
located in the vicinity of the hot star, within the region eclipsed 
by the giant. According to the modeling of the broad Ha wings 
(Sect. 3.2.2), the eclipsed emission region can be ascribed to the 
ionized wind from the hot star. 



4.1 .1 . Nebular lines 

Spectroscopic observations made after the major 1989-91 ac- 
tive phase at/around the 1994 eclipse showed that the nebular 
[Oiii] lines were subject to eclipse. Their fluxes faded by a fac- 
tor of ~2 with respect t o values meas ured before and after the 
eclipse (see Table 4 of ISkopal et al. 1 2001). In contrast, from 
2009 March, nebular [O m] line fluxes increased by a factor of 
~4.5 and kept their values at a constant level to our last 2009 ob- 
servation (2009/1 1/10), i.e. also throughout the eclipse (Fig. 4). 
This behavior implies that the hot star wind causes a fraction 
of the [Om] region. Its location is therefore a function of the 
mass-loss rate, which determines particle densities at a given 
distance from the wind source. Higher places the [O m] zone 
to longer distances from the hot star and vice versa. During the 
active phase, ~ 2 x 10"^ Moyr"', Vco - 2000 kms"' and 
critical densities for creation of nebular lines, ^ 10^ cm"^, 
place the [Om] zone at distances r > 210 Rq from the hot star, 
i.e. far beyond the eclipsed region. However, during quiescence, 
Mw ~ 1 X 10"^ Mqjv-^ and = 1600 kms"' (Table 7) cor- 
respond to the critical radius of ~50 Rq, which thus can be par- 
tially eclipsed by the giant with /?„ ~ 115 Rq. Finally, the [O iii] 
line profiles consist of more components (Fig. 3), which sug- 
gests that the hot star wind at farther distances from its origin 
was structured and not simply uniform in direction and density. 

4.2. The neutral zone 

We indicated the presence of a WSC of radiation in the com- 
posite spectrum photometrically by the 2009 eclipse profile, and 
spectroscopically by separating the observed spectrum into its 
individual components. Its source is characterized with effec- 
tive temperatures of ~5500 - 6500 K and effective radii of ~6 
-11 Rq, which correspond to luminosities between ~30 and 
100 Lq (Table 8). However, its shape cannot be spherical: (i) If 
this were a sphere, its radiation would not be capable of causing 
the observed nebular emission. On the other hand, the presence 
of the strong nebular component in the spectrum constrains the 
presence of a hot ionizing source in the system, (ii) If the radius 
of the eclipsed object, 7?e - 27 Rq (Sect. 3.1.2), were that of a 
sphere, its luminosity would be a factor of ~15 higher than that 
we observed at the 2009 bright stage. These arguments suggest 
that the source of the WSC, which is subject to the eclipse, has 
a form of a disk with the radius R^ - R^.. Assuming that its 
outer rim represents the warm pseudophotosphere (Sect. 5.3.5. 
of[Skopal 2005), then its luminosity 

InR^lH^crTt^ = AniRf^^JcrTt^. (11) 
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ACTIVE PHASE 




hot star 



Fig. 9. Sketch of the ionization structure of AX Per during the quiescent phase (left, pole-on view) and around the hot star during 
the active phase (right, edge-on view). During quiescence, the solid thick line represents the H'VH^ boundary in the steady-state 
approximation. Including the wind from the hot star and rotation of the binary (denoted by the arrow), a denser region evolves at 
th e preceding front of t he hot star, which is further twisted to the front of the giant against its motion (crosses, adapted from Fig. 5 
of lBisikalo etani2006h . The right panel shows a sketch of the ionization structure of the hot component during activity. The neutral 
zone has the form of a flared disk, whose optically thick rim represents the warm pseudophotosphere. The enhanced wind from the 
hot star is ionized by its radiation. The He+^ boundary is denoted by the gray dashed line (see Sect. 4). 



For Hy)IRy) = 0.3 (Sect. 3.2.2) and Z?^" = 11% (Table 8) we 
get Ry) - 20 Rq, which agrees well with R^ - 22.3 + 2 Rq, 
obtained from the first two contact times of the 2009 eclipse. 
Therefore we conclude that the warm stellar object, which de- 
veloped during the active phase, has the form of a flared disk, 
whose outer rim simulates the warm photosphere. Material lo- 
cated outside of the disk's rim is in the ionization shadow and 
thus represents the neutral zone (Fig. 9). 



4.3. Basic ionization structure 

The results we inferred above imply a significant difference 
in the ionization structure of the circumbinary environment of 
AX Per between its quiescent and active phase. The correspond- 
ing sketch is shown in Fig. 9. It can be described as follows. 

During the quiescent phase , the steady-state STB 
dSeaquist. Taylor. & Button] 1 19841) model is sufficient to 
explain the radio emission from symbiotic stars and the wave - 
like orbitally related variation in the optical LCs (ISkopalll200lh . 
However, to explain the shift of light minima to a position 
preceding the spectroscopic conjunction (Fig. 2, right panel), 
a more sophisticated model of the density distribution in the 
binary has to be considered. Hydrodynamical calculations of 
the structure of stellar winds in symbiotic stars that include 
effects of the orbital motion suggest a higher density region, 
twisted from t he preceding front of the hot star mo tion to that of 
the giant (e.g. lWalder]|1995l:lBisikalo et al. Il2006h . As a result, 
the light minima occur before the in ferior conjunction of the 
giant (see also Fig. 13 of lSkopanil99 8). In Fi g. 9 we plott ed a 
scheme of this region as calculated bv lBisikalo et al. I (12006) for 
a mass-loss rate from the giant and the hot star of 2 x 10"^ and 
1.2 X 10"^ Mq yr"', respectively. 

During the active phase, the neutral hydrogen zone has the 
form of a flared disk encompassing the hot star Its optically 
thick rim represents the warm pseudophotosphere. The ionized 
region fills in the remainder of the space above/below the neu- 
tral zone, which is ionized by the hot star It can be associated 
with the enhanced stellar wind from the hot star, whose emis- 
sion is subject to eclipse. The proposed ionization structure is 
consistent with the UV/near-IR model SED during active phases 
of symbiotic binaries with high orbital inclination (see Fig. 27 



in I Skopal I l2005h and that of the bipolar st ellar wind fro m the 
hot star used to model the broad Ha wings (ISkopal 112006) . The 
formation mechanism of the neutral disk-like zone is not under- 
stood yet. However, first results based on the theory of the wind 



ry . 

compression model dBiorkman & CassineUi 1993h support the 
findings of this paper dCarikova & Skopal 11201 lb . 



5. Conclusions 

We have investigated a new active phase of the eclipsing symbi- 
otic binary AX Per, which began during 2007. Our main conclu- 
sions may be formulated as follows. 

(i) AX Per entered a new active phase from ~2007.5 
(Sect. 3.1.1). After 10 orbital cycles (~ 18.6 years) we again 
measured a deep narrow minimum in the LC, which was caused 
by the eclipse of the hot component by its giant companion. We 
improved the ephemeris of eclipses (Eq. (1)) and found that their 
timing is identical (within uncertainties) with that of the inferior 
conjunction of the giant. We determined the radius of the giant as 
Rg = 115 + 2 Rq and that of the eclipsed object, 7?e = 27 ± 2 Rq 
(Sect. 3.1.2). 

(ii) The new active phase was connected with a signifi- 
cant enhancement of the hot star wind. From quiescence to 
activity, the mass-loss rate increased from ~ 9 x 10"^ to ~ 
2 X lO^^Moyr^', respectively (Sect. 3.2.2, Fig. 6, Table 7). 
The mean electron concentration within the eclipsed He^^ zone 
(^Heii ~ 5QRq above the hot star), ~ 10''^ cm"-', agrees well 
with that given independently by the mass loss rate (Sects. 3.2.1 
and 3.2.2). The hot star wind creates also a fraction of the [O iii] 
zone. Depending on the mass-loss rate, the nebular lines can be 
subject to partial eclipse (1994 eclipse) or not (2009 eclipse; 
Sect. 4.1.1). 

(iii) We revealed a significant change of the ionization struc- 
ture in the binary with respect to the quiescent phase (Sect. 4, 
Fig. 9). During the 2007-10 active phase we identified a warm 
stellar source located at the hot star equator (Sect. 3.2.3, Fig. 8). 
The source had the form of a flared disk, whose outer rim simu- 
lated the warm photosphere. It blocks out the ionizing radiation 
from the hot star, and consequently the material outside of the 
disk's rim is mostly in neutral form (Sect. 4.2). The ionized re- 
gion is located above/below the neutral zone and can be associ- 



11 



A. Skopal et al.: Formation of a disk structure in the symbiotic binary AX Per during its 2007-10 precursor-type activity 



ated with the enhanced hot star wind. The formation of the disk- 
like neutral zone and the increase of the nebular emission during 
the active phase was connected with the enhanced wind from the 
hot star (Tables 7 and 8). It is probable that this connection rep- 
resents a common origin of the warm pseudophotospheres that 
are observed during the active phases of symbiotic binaries. 
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Fig. 4. Variation in the selected line fluxes during the 2007-10 
active phase. Top panel shows the B-band LC to compare varia- 
tion in the continuum and to visualize better the time-interval of 
the eclipse (the light shadow band). The width of the dark and 
mild shadow band corresponds to the linear size of the eclipsed 
object and the He^^ zone, respectively. Dotted lines limit the ra- 
dius of the cool giant. Fluxes are dereddened and their values are 
given in Table 5. 
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Fig. 5. Radial velocities of main emission peaks measured for 
highly ionized elements (open symbols. Table 6) as a function of 
the photometric phase (1). Compare d are radial velocity curves 
for the M-type (iFekel et alTI l2000i solid line, filled circles) 
and A-type absorption lines (Mikol aiewska & Kenvon 1 1 19921 
dash ed hne, aster isks), and emission peaks of [Fevn]/16087 
line (lliiima 1119881 grey line, crosses). Systemic velocity of - 
117.44km s"' was subtracted. 
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Fig. 6. Example of the broad Ha wings from the 2007-10 active 
phase (solid line), compared with those observed during quies- 
cent phase (dashed line). Both spectra were taken at the same 
orbital phase, (f = 0.16. Local continua were subtracted from 
the spectra. The spectra were dereddened. The models (dotted 
lines) indicate a significant enhancement of the mass-loss rate 
during activity with respect to quiescence (Table 7, Sect. 3.2.2). 
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Fig, 7. Blue part of the dereddened spectra with well-developed 

warm pseudophotosfere as suggested by the model SED (Fig. 8). 
Features typical for the spectrum of F-G giants are marked by 
gray bands. 
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